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Abstract. We incorporate quasars into an analytic model 
to describe the reheating and reionization of the universe. 
In combination with a previous study of galaxies and 
Lyman-a clouds, we are able to provide a unified descrip- 
tion of structure formation, verified against a large variety 
of observations. We also take into account the clumping of 
the baryonic gas in addition to the presence of collapsed 
objects. 

We consider two cosmologies: a critical universe with 
a CDM power-spectrum and an open universe with Qq = 
0.3, A = 0. The derived quasar luminosity function agrees 
reasonably well with observations at z < 4.5 and with 
constraints over larger redshifts from the HDF. The radi- 
ation produced by these objects at z ~ 16 slowly reheats 
the universe which gets suddenly reionized at z r i = 6.8 
for the open universe [z r i = 5.6 for the critical density 
universe). The UV background radiation simultaneously 
increases sharply to reach a maximum of J 2 i = 0.18 at 
z = 2.6, but shows strong ionization edges until z < 1. 
The metallicity of the gas increases quickly at high z and 
is already larger than O.OIZq at z — 10. The QSO number 
counts and the helium opacity constrain the reionization 
redshift to be z r i ~ 6. We confirm that a population of 
faint quasars is needed in order to satisfy the observations. 
Due to the low reionization redshift, the damping of CMB 
fluctuations is quite small, but future observations (e.g. 
with the NGST) of the multiplicity functions of radiation 
sources and of the HI and Hell opacities will strongly con- 
strain scenarios in which reionization is due to QSOs. The 
reasonable agreement of our results with observations (for 
galaxies, quasars, Lyman-a clouds and reionization con- 
straints) suggests that such a model should be fairly real- 
istic. 



Key words: cosmology: large-scale structure of Uni- 
verse - galaxies: evolution - quasars: general - intergalactic 
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1. Introduction 

An important goal of cosmology is to describe the struc- 
ture formation processes which led to the wide variety of 
astrophysical objects we observe in the present universe, 
from Lyman-a clouds to galaxies and clusters. Several 
studies have shown that the usual hierarchical scenarios 
(like the standard CDM model) can provide predictions 
which agree reasonably well with observations for galaxies 
as well as for Lyman-a clouds. This corresponds to objects 
at z < 5. However, it is possible to constrain the earlier 
evolution of the universe by studying the reheating and 
reionization history implied by such models. Indeed, ob- 
servations show that the universe is highly photo-ionized 
by z = 5 and a large reionization redshift could imprint 
a signature on the CMB radiation. Moreover, future mis- 
sions such as NGST could for instance detect quasars at 
high redshifts z > 5. 

In this article, we present an analytic model for the 
reheating and reionization history of the universe, adopt- 
ing a CDM power spectrum in a critical density and in 
an open universe. Similar studies have been performed 
previously via numerical simulations (e.g. Gnedin & Os- 
triker 1997) and analytic approaches (e.g. Haiman & Loeb 
1997; Haiman & Loeb 1998) based on the Press-Schechter 
prescription (Press & Schechter 1974). However, previ- 
ous analytic models were often developed for this spe- 
cific purpose (i.e. they were not derived from a model 
already checked in detail against observations of galax- 
ies or Lyman-a clouds) and neglected the clumping of the 
gas (except for the presence of virialized objects used to 
count galaxies). Thus, the main motivations of our present 
study are to: 

- describe these early stages of structure formation 
through a self-consistent model which has already been 
applied to galaxies (Valageas & Schaeffer 1998) and to 
Lyman-a clouds (Valageas et al. 1999a). 

- take into account the broad range of density fluctua- 
tions within the IGM through our description of Lyman-a 
clouds. 
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- use this feature to constrain our model against several 
observations: notably the QSO number counts and the 
Gunn-Peterson test (for HI and Hell) . 

- develop a simple analytic model which can predict 
many properties of the universe (galaxy and quasar lumi- 
nosity functions, temperature and ionization state of the 
IGM, intensity and spectrum of the UV background radi- 
ation and fraction of matter within stars) and provide a 
complementary tool to numerical simulations. 

Consideration of the various objects involved in our 
work (beyond the just-virialized halos which are usually 
studied) is made possible because of a specific descrip- 
tion of the density field based on the assumption that the 
many-body correlation functions obey the scaling model 
detailed in Balian & Schaeffer (1989) and checked numeri- 
cally in Colombi et al.(1997). This allows one to define the 
various mass functions of interest, as described in Valageas 
& Schaeffer (1997; also in Valageas et al. 1999b), and to go 
beyond the scope of the usual Press-Schechter approxima- 
tion (Press & Schechter 1974). The main advantage of our 
approach is thus to provide a globally consistent picture of 
structure formation in the universe, within the framework 
of a hierarchical scenario. 

This article is organized as follows. In Sect. 2 we de- 
scribe our prescription for mass functions. Next, in Sect. 3 
we review our model for galaxy formation, described in 
more detail in Valageas & Schaeffer (1998) while in Sect. 4 
we deal with our prescription for quasars. In Sect. 5 we 
summarize the relevant aspects of our model for Lyman-a 
clouds (Valageas et al. 1999a). We describe the calculation 
of the evolution of the IGM properties (temperature, UV 
background radiation, ionization state) in Sect. 6 and in 
Sect. 7. Finally, in Sect. 8 and Sect. 9 we present our results 
for the case of an open universe and then for a critical 
universe. 

2. Multiplicity functions 

We first review the method we use to obtain the mass func- 
tions of various astrophysical objects, specifically galaxies 
and Lyman-a clouds. We consider objects of dark matter 
mass M to be defined by a density threshold A(M,z). 
This constraint depends on the class of astrophysical ob- 
jects one considers and it allows us for instance to distin- 
guish clusters from galaxies which correspond to higher 
density contrasts (see VS II). Lyman-a clouds are also 
formed by several populations of different objects which 
are not always defined by a constant density threshold 
(see Sect. 5). Note that such a goal is beyond the reach of 
the usual Press-Schechter prescription (Press & Schechter 
1974) which only deals with "just-collapsed halos" while 
we wish to describe simultaneously a wide variety of ob- 
jects. In any case, we attach to each halo a parameter x 
defined by: 

1 + A(M, z) 



where 

1{R) = 



d 3 n d 3 r 2 . . 4 3 

y2 £ ' Wlth V = 3 nR 



is the average of the two-body correlation function 
62 (ri , r2) over a spherical cell of radius R and provides the 
measure of the density fluctuations in such a cell. Then, 
we write the multiplicity function of these objects (defined 
by the constraint A(M, z)) as (see VS I): 



/, r x dM p , . . dx 
M M w x 



(2) 



where p is the mean density of the universe at redshift z, 
while the mass fraction in halos of mass between M and 
M + dM is: 



/, , x dM , . . dx 



(3) 



The scaling function H(x) depends only on the initial 
spectrum of the density fluctuations and must be obtained 
from numerical simulations. However, from theoretical ar- 
guments (see VS I and Balian & Schaeffer 1989) it is ex- 
pected to follow the asymptotic behaviour: 



x < 1 : H(x) cx x 



w-2 



x > 1 : H(x) oc x' 



0J3 — 1 p —x/x, 



with lo ~ 0.5, uj s ~ —3/2, x* ~ 10 to 20 and by definition 
it must satisfy 



f 

Jo 



x H(x) dx = 1 



(4) 



x{M,z) = 



£[R(M,z),z] 



(1) 



The correlation function £, that measures the non-linear 
fluctuations at scale R can be modelled in a way that 
accurately follows the numerical simulation. The mass 
functions obtained from (2) for various constraints A(M) 
were checked against the results of numerical simulations 
in Valageas et al. (1999b) in the case of a critical uni- 
verse with an initial power-spectrum which is a power-law: 
P{k) cx k n with n = 0,-1 and —2. This study showed 
that this model provides a reasonable approximation to 
the mass functions obtained in the simulations and that 
it works quite well for the two cases we shall need in the 
present article: i) a constant density threshold A ~ 178 
and ii) a constant radius constraint (or (1 + A) cx M). 
Moreover, the results of Valageas et al. (1999b) showed 
that H(x) is close to a similar scaling function h{x) ob- 
tained from the counts-in-cells statistics, as expected from 
theoretical considerations (see for instance VS I). 

It is clear that the model outlined above provides a 
unified description of various astrophysical objects which 
are obtained from the same non-linear density field. This 
is a great advantage of this approach since it ensures that 
we can model a wide variety of objects, from low density 
Lyman-a clouds to high density bright galaxies, in a fully 
consistent way. Then, we can study the interplay between 
these various structures as they develop progressively. 
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3. Galaxy formation 

In this paper we wish to study the reionization history 
of the universe. Since a large part of the ionizing radia- 
tion will be emitted by stars, we first need to devise a 
model for galaxy and star formation. We shall use a sim- 
plified version of the model described in detail in VS II 
and which was there compared with many observations. 
One can define galaxies by the requirement that two con- 
straints be satisfied by the underlying dark matter halo: 
1) a virilization condition A > A c (where A c {z) ~ 178 
is given by the spherical model and is constant for a criti- 
cal universe) and 2) a cooling constraint t coo i < tn which 
states that the gas must have been able to cool within a 
few Hubble times at formation. However, at high redshifts 
z > 1 the cooling constraint becomes irrelevant since any 
object which satisfies 1) also satisfies 2). Hence since we 
are mainly interested in large redshifts 2>lwe shall sim- 
ply define galaxies by the virialization condition A = A c . 
We also require that the virial temperature T of the halo 
be larger than the "cooling temperature" T coo i(z) at red- 
shift z. The latter corresponds to the smallest virializcd 
objects which can cool efficiently at redshift z, defined by 
the constraint: 



tcool — S iff 



(5) 



where s = 6 is a proportionality factor (one must have 
s > 1 since cooling is more efficient within the halo where 
the density is larger than on its boundary and cooling ac- 
celerates as baryons collapse). Here tu{z) is the age of the 
universe at redshift z while t coo i is the cooling time of a 
halo with density contrast A c (z), mass M, taking into ac- 
count both cooling (recombination, molecular cooling) and 
heating (by the background UV flux) processes. Since the 
physical properties of virialized halos with temperature T 
and density contrast A c (z) are different from the IGM, 
we let the chemical reactions (involving HI, HII, H~, H 2 , 
H^, Hel, Hell, HcIII and e~) evolve for a Hubble time tjy 
within this environment (defined by T and A c ) before we 
evaluate the cooling time t coo i. The main effect is that at 
large redshifts such clouds may produce enough molecular 
hydrogen to make molecular cooling efficient while with 
the use of the IGM abundances one would underestimate 
this contribution; see for instance Tegmark et al.(1997) for 
a detailed discussion. This will also appear clearly below 
in Fig. 4 where we compare the main contributions to cool- 
ing for both the IGM and these cooling halos. The virial 
temperature T coo i also defines the mass M coo i(z) and the 
radius R coo i{z) of the smallest objects which can cool and 
eventually form stars at redshift z. From the lower-bound 
T CO oi{z) and the virialization constraint A = A c (z), we 
obtain the mass function of galaxies at redshift z using 
(2). 

Next, we must attach a specific stellar content to these 
galactic halos. We shall again use the star formation model 
described in VS II. This involves 4 components: (1) short- 



lived stars which are recycled, (2) long-lived stars which 
are not recycled, (3) a central gaseous component which 
is deplenished by star formation and ejection by super- 
novae winds, and replenished by infall from (4) a diffuse 
gaseous component. The star-formation rate dM s /dt is 
proportional to the mass of central gas with a time-scale 
set by the dynamical time. The mass of gas ejected by 
supernovae is proportional to the star-formation rate and 
decreases for deep potential wells as 1 /T, in a fashion simi- 
lar to that adopted by Kauffmann et al.(1993). It was seen 
in VS II that for such a model a good approximation for 
the star- formation rate is: 



dM. 



= with 
dt To 



T * 1 + 



Tsn \ 
T ) 



(6) 



where M g is the total mass of gas, Td is the dynamical 
time and Tsn — 10 6 K describes the ejection of gas by 
supernovae and stellar winds: 



Tsn = 



2 e E S n Lim P Vsn 
3 k msN 



10 b K 



(7) 



Here e ~ 0.1 is the fraction of the energy E$n delivered by 
supernovae transmitted to the gas (Esn = 10 51 erg) while 
Vsn /msN — 0.005 Mq 1 is the number of supernovae per 
solar mass of stars formed. Note that for halos defined 
by a constant density threshold A c ~ 178 we have Td ~ 
tff(z). Although (6) was obtained for small galaxies with 
T <C Tsn (which is the range we are mainly interested 
in) it also provides a reasonable approximation for large 
galaxies T > Tsn- In the case f2 = 1 we obtain in our 
model for a galaxy similar to the Milky Way (i.e. with 
a circular velocity V c = 220 km/s): t ~ 7 10 9 years and 
dM s /dt ~ 5M /year (see VS II). This star formation rate 
is consistent with observations (McKee 1989). Then the 
mass of gas at time t within the galaxy is given by: 



M g = M, 



go e 



-t/r 



(8) 



where M g o = Mb is the initial mass of baryons which we 
take to be proportional to the dark matter mass M: 



M h = ^M 



n 

From this model, the star-formation rate per Mpc 3 is: 



(9) 



dp s 
dt 



fi& p(z) 



P_ 

Pd 



T 



SN 

T 



(10) 



x e 



' x H(x) — 

X 



where p/Pd = 0.7 is a parameter of order unity which 
enters the definition of the dynamical time Td. The signifi- 
cance of each term in this expression is clear and the tem- 
perature dependence simply states that the average star 
formation efficiency of small galaxies is small as the gas 
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is easily expelled by supernovae. Note that in the origi- 
nal model described in VS II for bright galaxies at low 
redshifts, the star formation rate declines since most of 
the gas has already been consumed. This does not ap- 
pear in (10) because we defined all galaxies by A = A c 
while at low z for large T the cooling constraint implies 
that (1 + A) cx M (i.e. R is constant) which decreases the 
galactic dynamical time and increases the ratio t/r which 
enters (8). 

To derive the radiation emitted by galaxies, we do not 
need their global star formation rate but their stellar con- 
tent. However, as shown in VS II, the mass in the form of 
short-lived stars (i.e. with a life-time T s h small compared 
to tn) of mass to to to + dm is given by: 



dM sh = dri — M g = drj T sh ( 

To \ dt 



(11) 



where drj = mcf){m)dm is the fraction of mass which goes 
into such stars for each unit mass of stars which are 
formed. This depends on the initial stellar mass function 
(IMF) 4>(m). Since the stellar radiation output at high 
energy (y > 13.6 eV) is dominated by the most massive 
stars, the relation (11) will be sufficient for our purposes. 
Next, if we assume that stars radiate as blackbodies with 
an effective temperature T e ff cx L 013 and we use the 
mean scalings r s h cx m/L and L cx to 3,3 we obtain the 
energy output of such galaxies: 



d 2 E \ 
dvdt) 



dM s 
~~dt 



lyr 



(12) 



with 



L vs (v) = — / 



m(j)(m)dm 



2-Khv 4 



aT A c 2( e hv/kT _ X ) 



From the radiation emitted by individual galaxies we now 
wish to estimate the energy received by a random point 
in the IGM. We shall write the source term S us due to 
stellar radiation for the background UV flux J„, see (30), 
as the following average: 



(13) 



where r\g{x)dx j x is the mass function of galaxies, obtained 
from (2) as described previously, while r s is a mean opacity 
which takes into account the fact that the radiation emit- 
ted by galaxies can be absorbed by the IGM and Lyman-a 
clouds. We shall come back to this term later. Thus, we 
get in this way a simple model for the stellar radiative 
output from our more detailed description of galaxy for- 
mation. The reader is referred to VS II for a more precise 
account of the details and predictions of our galaxy forma- 
tion model. Note that our prescription is consistent with 
such observations as the Tully-Fisher relation and the B- 
band luminosity function. 



4. Quasar radiative output 

In addition to galaxies we also need to describe the radi- 
ation emitted by quasars which provide a non-negligible 
contribution to the background radiation field, especially 
at the high frequencies v > 24.6 eV which are relevant for 
helium ionization. We shall again follow the formalism of 
VS II to obtain the quasar luminosity function, in a fash- 
ion similar to Efstathiou & Rees (1988) and Nusser & Silk 
(1993). We assume that the quasar mass Mq is propor- 
tional to the mass of gas M gc available in the inner parts 
of the galaxy: Mq = F M gc . Note that for galaxies which 
have not yet converted most of their gas into stars (i.e. all 
galaxies except those with T > Tsn at z < 1) this also 
implies Mq ~ F M s where M s is the stellar mass. Indeed, 
for tn < To (where tn is the age of the universe) we have: 



M s ~ t H /T Mg 



(14) 



by definition of tq, see (6), while the mass M gc of cold 
central gas M gc satisfies: 



T. 



SN 



M„ 



(15) 



The factor (1 + Tsn/T) translates the fact that in our 
model supernovae eject part of the star-forming gas out of 
the galactic center into the larger dark matter halo (VS 
II). Hence we get M gc ~ M s . Of course, at late times for 
bright galaxies when most of the gas has been consumed 
we have M gc <C M s . Then the mass of gas available to 
feed the quasar declines with time. This leads to a high 
luminosity cut-off at low z for the quasar luminosity func- 
tion since in this regime very massive galaxies have less 
gas than smaller ones which underwent less efficient star 
formation (see VS II and Sect.8.1). We shall use F = 0.01 
for fl = 1 and F = 0.006 for fl = 0.3. Note that observa- 
tions (Magorrian et al.1998) find Mq ~ 0.006 M s in large 
galaxies. Next we write the bolometric luminosity Lq of 
the quasar as: 



£ Mq C 2 
tQ 



(16) 



where e = 0.1 is the quasar radiative efficiency (fraction of 
central rest mass energy converted into radiation) and tQ 
is the quasar life-time. Since we shall assume that quasars 
radiate at the Eddington limit we have: tQ = 4.4 e 10 8 yr. 
Thus, the quasar luminosity attached to a galaxy of dark 
matter mass M, virial temperature T, is: 



eF fl b 

tQ Qq 



Tsn \ 
T ) 



Mc 2 



(17) 



As seen above in (15), the temperature term comes from 
the fact that in our galactic model, small galaxies (T <C 
Tsn) are strongly influenced by supernova feedback which 
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expells part of their baryonic content from the inner re- 
gions. Note however that this term does not enter the 
relation (quasar mass) - (stellar mass) as it cancels out on 
both sides. Next we obtain the quasar multiplicity func- 
tion from the galaxy mass function as: 



, , , n dM Q , , .> dM 



tM 



The factor Aq < 1 (we use Aq = 0.06) is the fraction 
of galactic halos which actually harbour a quasar while 
t m is the evolution time-scale of galactic halos of mass M 
defined by: 



(19) 



Since the quasar life-time ig ~ 10 8 yr is quite short, this 
reduces to tjq (Mq ) <1Mq / Mq = Xq t Q dr) g /dt dM/M. 
Together with (17) the relation (18) provides the quasar 
luminosity function. Thus, we only have two parameters: 
(e F/to) (which only depends on F, constrained by the 
observed (quasar mass)/(stellar mass) ratio, for quasars 
shining at the Eddington luminosity) which enters the 
mass-luminosity relation, and (Aq £q) which appears as 
a simple normalization factor in the luminosity function. 
Hence a larger fraction of quasars Aq together with a 
smaller life-time tQ would give the same results, so that 
we could also choose Aq = 1. In a fashion similar to what 
we did for galaxies we can now derive the quasar radia- 
tive output. We first write the radiation emitted by an 
individual quasar as: 



d 2 E 
dvdt 



Lq 

VB 



L B 
Lboi 



(20) 



where (Lb/L^) = 0.094 is the conversion factor from 
bolometric luminosity to B-band luminosity (Lb = 
vbL v (vb) at vb — 2.8 eV), taken from Elvis et al.(1994), 
while a = 1.5 is the local slope of the quasar spectrum. 
Then, the source term S u q for the background radiation 
due to quasars is: 



dx 



fd 2 E\ 



(21) 



where again tq is a mean opacity which we shall describe 
later. 



5. Lyman-a clouds 

The description of gravitational clustering used in this 
article allows one to build a model for Lyman-a clouds 
(Valageas et al. 1999a). We shall take advantage of this 
possibility to include these objects in the present study. 
Indeed, although at high rcdshifts they do not contribute 
significantly to the total opacity (which comes mainly 
from the uniform component of the IGM) since only a 



small fraction of baryonic matter has been allowed to form 
bound objects, at redshifts close to the reionization epoch 
they already provide a non-negligible opacity. We identify 
Lyman-a absorbers as three different classes of objects, 
which we shall briefly describe below. 



(18) 

5.1. Lyman-a forest 



We assume that after reionization the gas within low- 
density halos is reheated by the UV flux to a temperature 
Tl v = 3 10 4 K. Hence in such shallow potential wells, 
baryonic density fluctuations are erased over scales RdLy 
defined as in (22) but with the temperature Tl v . This 
builds our first class of objects defined by their radius 
RdLy(z) and virial temperatures T < Ti y . The multiplic- 
ity function of these mass condensations is again obtained 
from (2). The fraction of neutral hydrogen at low z is eval- 
uated by assuming photo-ionization equilibrium. At high 
z prior to reionization, when the UV flux is very small and 
cannot heat the gas, we simply take T^ y = Tiqm while 
the fraction of neutral hydrogen is unity. Since the bary- 
onic density is roughly uniform within these objects (by 
definition) we consider that each halo produces one spe- 
cific mean column density on any intersecting line-of-sight 
(we neglect the small dependence on the impact parameter 
due to geometry). At low z this population can be identi- 
fied with the Lyman-a forest. Note that, as explained in 
details in Valageas et al. (1999a), our approach is also valid 
for clouds which are not spherical objects of radius RdLy 
but filaments of thickness RdLy and length L 3> RdLy 
This is due to the growth of the density fluctuations on 
smaller scales (along with £) and to the direction jumps 
of filamentary structures. 

Here we note that models for the Lyman-a forest are 
often classified in two categories: 1) mini-halo models and 
2) IGM density fluctuations. In case 1), one considers that 
Lyman-a absorbers are discrete clouds formed by bound 
collapsed objects (or halos confined by the IGM pressure) 
which occupy a small fraction of the volume. On the other 
hand, in case 2) (which is currently favored) one assumes 
that absorption comes from a continuous medium (the 
IGM) with relatively small density fluctuations. Although 
in our model we identify distinct patches of matter (of 
size RdLy) as in 1), the underlying picture corresponds 
to the case 2). Indeed, as we consider regions with an 
"overdensity" (1 + A) from <~ 20 down to (1 + A)jgm, 
defined below in (24), which can be as low as 10~ 3 , see 
Fig. 13, we take into account all the volume of the universe. 
Hence our Lyman-a forest absorbers are made of a broad 
range of density fluctuations within the IGM which fill 
all the space between galactic halos (which we describe 
below as they form Lyman-limit and damped systems and 
only occupy a negligible fraction of the volume, as seen in 
Fig. 12). Note that this would not be possible if we were to 
consider density fluctuations defined by a constant density 
threshold (1 + A) t h > 1 since this would imply that we 
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probe at most a fraction 1/(1 + A) t h of the volume of 
the universe. We identify the lowest density regions (i.e. 
with a density contrast A;gm), which are also the most 
numerous and fill most of the volume, with the IGM. A 
patch of matter with this density would only make up a 
column density Nhi ~ 10 6 cm~ 2 on a scale RdLy at z = 0. 

5.2. Lyman-limit systems 

Potential wells with a large virial temperature T > T^ y 
do not see their baryonic density profile smoothed out 
and they also retain their individuality. Thus, we define a 
second class of objects identified to the ionized outer shells 
of virialized halos, characterized by their density contrast 
A c and satisfying T > T^ y . The deepest of these potential 
wells (such that T > T coo i) corresponds to the galactic 
halo described in Sect. 3. We assume that the mean density 
profile is a power-law p oc r~ 7 (with 7 = 1.8) so that each 
object can now produce a broad range of absorption lines, 
as a function of the impact parameter of the line-of-sight. 
This population can be identified with the Lyman-limit 
systems. 



5.3. Damped systems 

The deep cores of the virialized halos described above 
which are not ionized because of self-shielding (at low z) 
form our third population of objects. One halo can again 
produce various absorption lines for different impact pa- 
rameters. At high z, prior to reionization, halos are en- 
tirely neutral so that the previous contribution of ionized 
shells disappears and we only have two classes of objects: 
these neutral virialized halos and the "forest" objects. 



6. Evolution of the IGM 

We now turn to the IGM itself. We model the universe at 
a given redshift z as a uniform medium (the IGM), char- 
acterized by a density contrast Aiqm, a gas temperature 
Tigm and a background radiation field J v , which contains 
some mass condensations recognized as individual objects 
identified with galaxies or Lyman-a clouds as described 
above. 

Since the gas in the IGM has non-zero temperature 
Tigm, baryonic density fluctuations are erased over scales 
of order Rd{z) within shallow potential wells with a virial 
temperature T vir < Tigm or within "voids" , with: 



Note that the damping scale Rd is different from the Jeans 
scale: 



1 / jkT IGM 

2 H V pm p 



(22) 



where C s is the sound speed, tn the age of the universe, 
m p the proton mass and 7 <~ 5/3. Indeed, the pressure 
dominates over gravitation for objects such that T vir < T. 



Rj{z) 



4irGpm p p D M 



(23) 



Both scales are equal (up to a normalization factor of or- 
der unity) if the dark matter density is equal to the mean 
universe density: pdm = P- However, we shall consider 
underdense regions where (1 + A) can be as low as 10~ 3 , 
see Fig. 13 below. Indeed, as an increasingly large propor- 
tion of the matter content of the universe gets embedded 
within collapsed objects as time goes on the density of 
the IGM (the volume between these mass condensations) 
becomes much smaller than the mean universe density. In 
this case where pdm < p we have Rd < Rj. We use Rd 
because of the finite age of the universe: the medium can- 
not be homogenized over scales larger than those reached 
by acoustic waves over the time tn (the scale Rj corre- 
sponds to the limit of large times). Note that for Lyman- 
a clouds we also use Rd as the characteristic scale, with 
Tl v = 3 10 4 K, since we consider regions with very low or 
moderate densities (1 + A) < 45, see Sect. 5 and Valageas 
et al. (1999a). Then, the density contrast of the IGM is 
given by: 



(1 + A.) igm = Min 1 , £(R d ) 



0/(1- 



(24) 



This simply states that at high z (when £(R d ) <C 1) 
we have Pjgm = ~P (he. the universe is almost exactly 
a uniform medium on scale R d ) while at low z we have 
Pigm < P since most of the matter is now within over- 
dense bound collapsed objects (clusters, filaments etc.) 
while most of the volume (which we call the IGM) is 
formed by underdense regions. 

Since the mean density of the universe is < p >= ~p we 
define a baryonic clumping factor C& =< p\ > / < pb > 2 
by: 



/dr 
(1 + A) x 2 H(x) — 

~ (1 + A) 2 GM + F Lv < 1 + A > Ly +F mr (l + A c ) 



(25) 



where we used the fact that the volume fraction Figm,voI 
occupied by the IGM is very close to unity. Here Fl v and 
F vir are the fractions of mass formed by Lyman-a forest 
clouds (with a density contrast lower than A c ) and by 
virialized objects. Note that C& somewhat underestimates 
the actual clumping of the gas since we did not take into 
account the collapse of baryons due to cooling nor the 
slope of the density profile within virialized halos. How- 
ever, these latter characteristics are included in our model 
for Lyman-a clouds. We also define the mean density due 
to objects which do not cool as: 



< 1 + A>„= (1 + A) /GM + / x 2 H(x)— (26) 

Jo 



x 
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Before reionization this corresponds to the density field 
of neutral hydrogen since galactic halos (i.e. massive po- 
tential wells with x > x coo i which can cool) ionize most of 
their gas because of the radiation emitted by their stars or 
their central quasar. We obtain the mean square density 
in a similar fashion: 



< (1 + A) 2 >„= (1 + A) 



IGM 



Xcoal J 



and the corresponding clumping factor is simply: 



< (1+A) 2 >r 

< 1 + A >l 



(27) 



The quantities < 1 + A >„ and C n characterize the den- 
sity fluctuations of neutral hydrogen within the IGM. Note 
that most of the volume is occupied by regions which sat- 
isfy (1 + A) ~ (1 + A) IGM - 

The gas which is within the IGM is heated by the UV 
background radiation while it cools due to the expansion 
of the universe and to various radiative cooling processes. 
Note that we neglect here the possible heating of the IGM 
by supernovae. However supernova feedback is included 
in our model for galaxy formation: we simply assume it 
only affects the immediate neighbourhood of these galax- 
ies (see also McLow & Ferrara 1998). Thus, we write for 
the evolution of the temperature of the IGM: 



dTi G M „ a Tiqm , Tjcm 

-J = -2 - llGM ~ — + 



(28) 



where a(t) is the scale factor (which enters the term de- 
scribing adiabatic cooling due to the expansion) . The heat- 
ing time-scale theat is given by: 



heat 3/2n b kT IGM 



J2 J nwWly - v j) J »^ 



(29) 



where j = (HI,HeI,HeII), fj is the ionization threshold 
of the corresponding species, rij its number density in 
the IGM and nt, the baryon number density. The cool- 
ing time-scale t coo i describes collisional excitation, colli- 
sional ionization, recombination, molecular hydrogen cool- 
ing, brcmsstrahlung and Compton cooling or heating (e.g. 
Anninos et al.1997). Next, we can write the evolution 
equation for the background radiation field J u : 

^T = -3- Ju + -v^r- - KJ V +S vs + S vQ (30) 
at a a ov 

The first two terms on the r.h.s. describe the effects of 
the expansion of the universe, while the last two terms 
represent the radiation emitted by stars and quasars which 
we obtained previously. The absorption coefficient k v is 
written as: 



IMpc 



ijv,IGM + T l,NHI + T l,NHeI + T v,NHeIl) (31) 



where t 1 . IGM is the opacity at frequency v of the IGM over 
a physical length of 1 Mpc, while T v ^ j corresponds to the 
contribution by "Lyman-a" clouds (i.e. discrete mass con- 
densations as opposed to the uniform component which 
forms the IGM). Thus we have: 



T l.IGM = \ 'YjVA l ') n j 1M P C 



(32) 



Note that in this study we consider the medium as purely 
absorbing and we neglected the reprocessing of ionizing 
photons. From the evolution of the IGM temperature 
and the background radiation field we can also follow the 
chemistry of the gas within this uniform component. More 
precisely we consider the following species: HI, HII, H~, 
H2, H^~, Hel, Hell, Helll and e~ (see for instance Abel 
et al.1997 for rate coefficients). Thus we obtain the reion- 
ization history of hydrogen and helium together with the 
spectral shape of the background radiation J v . 

7. Opacity 

In the previous calculations, (13) and (21), where we de- 
scribed the source terms for the radiation field within the 
IGM we introduced opacity factors to model the absorp- 
tion of the radiation emitted by quasars and stars by the 
IGM and Lyman-a clouds. We shall deal with these terms 
in this section. 

We consider that each source (galaxy or quasar) active 
at a given redshift z ionizes its surroundings over a radius 
Ri given by: 



Ri 



HII 



4naC n n 2 



H 



-aC n n H t\ dN. t 
' dt 



1/3 



(33) 



where a(3 10 4 K) is the recombination rate (within the 
ionized bubbles), uh the mean number density of hydro- 
gen obtained from (26) , C„ the clumping factor from (27) , 
(dNj/dt) the emission rate of ionizing photons from the 
source and t its age. 

We have t < tn so we neglected here the influence of 
the expansion of the universe and the timc-dcpcndcncc 
of the source luminosity over its age. We take t = tn for 
galaxies and t = tQ for quasars. Although this procedure is 
consistent with our prescription for the quasar luminosity 
function (we assumed quasars to shine at the Eddington 
luminosity on the time-scale tQ <C tn and then to fade) 
we somewhat overestimate the radiative output of galax- 
ies since the galaxy luminosity function decreases with z 
over the time-scale tn- However, the relation (33) should 
still provide a correct estimate of the magnitude of this ef- 
fect. Note that Ri is smaller than the usual "Stromgren" 
radius which corresponds to the limit t — > 00 in (33). In- 
deed, the exponential term in (33) can also be written as 
exp(— t/t rec ) which shows that at redshifts z <C 100 where 
the recombination time is larger than the age of the source 
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(which is smaller than tn) the ionization front is smaller 
than the Stromgren radius. This effect was also described 
by Shapiro & Giroux (1987). In addition, these authors 
took into account the expansion of the universe but as- 
sumed a fixed number of sources. Here since we consider 
sources with a life-time t < tn we neglect the influence 
of the expansion of the universe but we take into account 
the increasing number of galaxies and quasars. 

Next, we obtain the volume fraction Qhii 3 q (i-e. the 
filling factor) occupied by such ionized bubbles around 
galaxies or quasars as: 



QhII 3 , q = J Vbubble 3 , Q (x) 



dx 4ir 3 



(34) 



For bubbles ionized by stellar radiation we have 
Vbubbie s ( x ) = Vg( x ) where r) g {x)dx/x is the mass function 
of galaxies while for quasars we write: 



Vbubbie Q {x) = Xq %(x) Min 



l,Max 



(tQ tree \ 
thl ' tM ) 



(35) 



where t rec is the recombination time within the ionized 
bubbles. This differs from the quasar mass function (18) 
through the term t rec /tM because a region remains ion- 
ized over a time-scale t rec which may be longer than the 
quasar life-time tq . Note that our general procedure only 
provides an upper bound to the actual efficiency of radia- 
tive processes since we did not include absorption within 
the host galactic halo itself. We do not integrate Qhii SiQ 
over time since this is already done in (33) and the sharp 
rise with time of the luminosity functions (before reion- 
ization) ensures that the radiative output is dominated 
by recent epochs. Moreover, since we have t rec < tn, as 
shown below in Fig. 8 (curve t rec ^ubbie) , ionized bubbles do 
not survive more than a Hubble time unless new sources 
(galaxies or quasars) appear. The filling factor within the 
IGM is written as the sum of the contributions from galax- 
ies and quasars: 



Qhiijgm — Qhii s + QhIIq 



(36) 



Of course the previous considerations only apply to high 
redshifts prior to reionization when the universe is almost 
completely neutral. At reionization these bubbles overlap 
and the background UV flux gets suddenly very large as 
absorption drops. At later times the whole universe is ion- 
ized so there are no more discrete bubbles (and formally 
Qhii — I). We also define in a similar fashion the filling 
factors QhcII and QhcIH which describe regions around 
quasars where helium is singly or doubly ionized. Since the 
stellar radiation shows an exponential decrease at high 
frequencies quasars are the only relevant source for this 
process. The filling factor Qhiijgm obtained above will 
be used to obtain the IGM opacity. 

However, as we explained previously we also consider 
the universe to contain numerous clouds which contribute 
to the opacity seen by the radiation field. We shall first 



consider that these clouds are ionized (or more exactly 
that their number density drops significantly) within the 
radius Ri, c i = Ri defined in (33) from the quasar. In 
other words, most of the opacity comes from clouds lo- 
cated deeply within the IGM where the background ra- 
diation J v is very small (before reionization) since close 
to quasars (within Ri >c i) the local radiation suddenly 
gets much higher. Note that the distribution of Lyman-a 
clouds we calculate is indeed obtained from the IGM back- 
ground radiation, which calls for the cutoff Ri iC i- However, 
at low z after reionization the "sphere of influence" of a 
quasar is no longer given by Ri, c i (since the whole medium 
is ionized). Instead, we define the radius Rj yC i by: 



/ ciiV, 



V dt 



10 



■h 



(37) 



where h is Planck constant and J21 is a measure of the 
background radiation within the IGM in the ionizing part 
of the spectrum: 



J. 



21 



JJv<THl(v) % 



(38) 



Thus, the "sphere of influence" of a quasar is defined by 
the region of space around the source where the radiation 
emitted by this quasar is significantly larger than the back- 
ground radiation (at high z when the medium is neutral 
this corresponds to Ri^i, while at low z when the IGM is 
ionized this is given by Rj c i)- Thus, in practice we shall 
simply use R c i = Min[i?i. c ;, R,j, c i] to obtain the volume 
fraction Qhii, cloud where the number density of Lyman-a 
clouds is significantly lower than within the IGM: 



„ f dx At: 

QHILcloud — / Vbubble Q (X) — "g"-^, 



(39) 



As we shall see from the numerical results, at low red- 
shifts z < z r i when the universe is reionized the opacity 
is very low (the UV flux is large) so that absorption plays 
no role for the evolution of J„. Thus, in practice the ra- 
dius Rj. c i is irrelevant. It only provides some information 
on the properties of the universe but it does not influ- 
ence the behaviour of the latter. Thus, while Qhiijgm 
increases with time until it reaches unity as the universe 
gets reionized, Qnii.cioud will first grow before reioniza- 
tion as the volume occupied by the ionized bubbles in- 
creases and then decrease at z <C z r i because the quasar 
luminosity function drops at low z. Since a fraction of vol- 
ume Q translates into the same fraction Q along a random 
line of sight (neglecting correlations in the distributions of 
sources) we write the opacity T v {r) seen from a point in 
the IGM to a source located at the distance r as: 



T A r ) = X! ( T ijV 3 Qj, cloud + rl IGM . QjJGM^j 
1 



IMpc 



(40) 



where Qhi = 1 — Qhii is the neutral hydrogen filling fac- 
tor (Q H ei = 1-QHeii-QHein) and 7-1 N corresponds to 
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discrete clouds while IGM . describes the IGM contribu- 
tion. The typical distance l g (x) between galactic sources 
characterized by their parameter x, density contrast A 
and radius R, is given by their number density, see (2), 



i fl (x)~JZ(l + A) 



1/3 



l H{x)} 



-1/3 



(41) 



where we did not take into account correlations. Since only 
a fraction Aq < 1 of galaxies host quasars we have for the 
mean distance between bubbles ionized by quasars: 



Iq(x) = A Q Min 



l,Max I ^ ' r " 



in tn 



-1/3 



l g (x) (42) 



as in (35). Here t rec is the recombination time within the 
ionized bubbles, see Fig. 8. Since at low redshift t rec ~ t H 

— 1/3 

we usually have Iq ~ Aq l g . Next, we define an effective 
opacity r e f f over the region of size I and volume V by: 

d 3 r 



= -r(r) 



V 



(43) 



where r(r) is given by (40). Then we use for the opacities 
which enter the source terms (13) and (21) a simple pre- 
scription which recovers the asymptotic regimes r e / / — ► 
and T e ff — > oo of (43): 



T ',Q ( x ) = hi ( 1 + -T u (l gtQ ) 



1 



T A l g,Q)" 



(44) 



8. Numerical results: open universe 



We can now use the model we described in the previous 
sections to obtain the reionization history of the universe, 
as well as many other properties such as the population of 
quasars, galaxies or Lyman-a clouds, for various cosmolo- 
gies. We shall first consider the case of an open universe 
VLq = 0.3, A = 0, with a CDM power-spectrum (Davis et 
al.1985), normalized to eg = 0.77. We choose a baryonic 
density parameter fi& = 0.03 and H = 60 km/s/Mpc. 
We use the scaling function h(x) obtained by Bouchct 
et al.(1991) as explained in VS II. Our model is consis- 
tent with the studies presented in VS II and Valageas et 
al. (1999a), so that those papers are part of the same uni- 
fied model and describe in more details our predictions for 
galaxies and Lyman-a clouds at z < 5. 

8.1. Quasar luminosity function 

Although our model was already checked in previous stud- 
ies for galaxies and Lyman-a clouds as explained above, 
our description for quasars was not compared to observa- 
tions in great detail (although a first check was performed 
in VS II). Thus, we first compare in this section our predic- 
tions for the quasar luminosity function to observational 
data, as shown in Fig.l. 

We can see that our model is consistent with observa- 
tions. At low redshifts the number of quasars we predict 




-5 



O. 1 6<z<0.4 
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1 .9<z<3.2 



I I I — I I I I I I — I — H4— I — I I I 



S.5<z<3.5 




3.5<z<4.5 



log(L B /L s ) 

Fig. 1. The evolution with redshift of the B-band quasar lumi- 
nosity function in comoving Mpc~ 3 . The data points are from 
Pei (1995). 



does not decline as fast as the data, however we get a 
significant decrease which is already an improvement over 
the results of Efstathiou & Rees (1988) for instance. We 
can note that Haehnelt & Rees (1993) managed to ob- 
tain a good fit to the observed decline at low z but they 
had to introduce an ad-hoc redshift and circular velocity 
dependence for the black hole formation efficiency. Since 
our model appears to work reasonably well we prefer not 
to introduce additional parameters. Moreover, as we no- 
ticed earlier our ratio (black hole mass) /(stellar mass) is 
consistent with observations (F — 0.006) while the quasar 
life-time we use tg = 0.44 10 s yrs agrees with theoret- 
ical expectations. The high-luminosity cutoff, which ap- 
pears at z < 3.5, comes from the fact that in our model 
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very massive and bright galaxies have consumed most of 
their gas. Thus, the maximum quasar luminosity starts 
decreasing with time at low z because of fuel exhaustion. 
We note that Haiman & Loeb (1998) obtained similar 
results at z > 2 although they used a very small time- 
scale tQ ~ 6.6 10 5 yrs (in our case this problem is partly 
solved by the introduction of the parameter Aq < 1 which 
states that only a small fraction of galaxies actually host 
a black hole) . However, they note that the number density 
of bright quasars they get increases until z — 0. 

In a recent paper, Haiman et al.(1998) point out that 
the lack of quasar detection down to magnitude V = 29 
in the HDF strongly constrains the models of quasar for- 
mation, which tend to predict more than 4 objects (which 
is still marginally consistent). In particular, they find that 
one needs to introduce a lower cutoff for the possible mass 
of quasars (shallow potential wells with a circular velocity 
lower than 50 km/s are not allowed to form back holes) 
or a mass-dependent black-hole formation efficiency. We 
show in Fig. 2 the predictions of our model. 




2 6 28 
magnitude 



Fig. 2. The quasar cumulative V-band number counts. The 
dashed line shows the counts of quasars with magnitude 
brighter than V located at the redshifts 3.5 < z < 4.5 while 
the solid line corresponds to 3.5 < z < z ri . 



The solid line shows the number N(< V) of quasars 
with a magnitude lower than V located at a redshift larger 
than 3.5 up to the reionization redshift z r i = 6.8. The 
large opacity beyond this redshift prevents detection of 
higher z objects. We also take into account the opacity 
due to Lyman- a clouds at lower z. We can see that our 
model is marginally consistent with the constraints from 
the HDF since it predicts 4 detections up to V = 29. We 
note that our model automatically includes photoioniza- 
tion feedback (threshold T coo i) and a virial temperature 
dependence in the relation (black hole mass) - (dark mat- 
ter halo mass). However, the "cooling temperature" T coo i 
is too low to have a significant effect on the number counts. 
Of course, we see that at bright magnitudes most of the 
counts come from low-redshift quasars (z < 4.5). Thus, 



the QSO number counts strongly constrain our model 
since in order to obtain a reionization history consistent 
with observations (namely the HI and Hell Gunn-Peterson 
tests and the low-redshift amplitude of the UV back- 
ground radiation field) we need a relatively large quasar 
multiplicity function. However, one might weaken these 
constraints by using an ad-hoc QSO luminosity function 
with many faint objects (Lb < 5 10 9 Lq). 

8.2. Reheating of the IGM 

As we explained previously the radiation emitted by galax- 
ies and quasars will reheat and reionize the universe, fol- 
lowing (28) and (30). We start our calculations at zi — 200 
with the initial conditions used by Abel et al.(1998), see 
also Peebles (1993). In particular: 



TwM(zi) = Min 



135 



1 



100 



K , 2.73 (1 + Zi ) K 



2.4 lo -4 n 



1/2 



0.05 



2 10 



-20 



(i-r)fl g 



3/2 



(1 + Zi) 



5.1 



and we use a helium mass fraction Y = 0.26. We present 
in Fig. 3 the redshift evolution of the IGM temperature 
Tig m , as well as the temperature T coo i which defines the 
smallest virialized objects which can cool at redshift z, see 
(5). 



« 4 





III 


1 1 1 | 1 1 1 


1 


















- _ _ x xv ~^ i 


















IGM 


\\ \ 

\ s 








\ \ 










\ \ 








,,, 


\ \ 







log(l +z) 

Fig. 3. The redshift evolution of the IGM temperature Tigm 
(solid curve) . We also show the virial temperature T coo i of the 
smallest virialized halos which can cool at redshift z (dashed 
curve) while T m is a mass-averaged temperature (dot-dashed 
curve) . 



At high z the IGM temperature decreases with time 
due to the adiabatic expansion of the universe. Next, for 
z < 24 (log(l + z) < 1.4) the medium starts being slowly 
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reheated by the radiative output of stars and quasars until 
it reaches at z ~ 9 a maximum temperature T max ~ 3 10 4 
K where collisional excitation cooling is so efficient that 
the IGM temperature cannot increase significantly any 
more. As we shall see later, this phase occurs before the 
medium is reionized, as was also noticed by Gnedin & 
Ostriker (1997) using a numerical simulation. There is a 
small increase at z ~ 7 (log(l + z) ~ 0.9) when the uni- 
verse is fully reionized and the UV background radiation 
shows a sharp rise. However, because cooling is very ef- 
ficient, the dramatic increase in J v only leads to a small 
change in Tigm- Eventually at low redshifts the temper- 
ature starts decreasing again due to the expansion of the 
universe as the heating time-scale becomes larger than the 
Hubble time tu . The temperature T coo \ which defines the 
smallest objects which can cool at redshift z increases with 
time because the decline of the number density of the var- 
ious species, due to the expansion of the universe, makes 
cooling less and less efficient. Indeed, the cooling rate (in 
erg cm -3 s~ 4 ) associated with a given process involving 
the species i and j can usually be written as kij(T)nirij, 
which leads to a cooling time-scale: 

t CO oi,n = * » _1 - (1 + tn (45) 

where we have neglected the temperature dependence. 
Thus, the ratio (cooling time) /(Hubble time) increases as 
time goes on (at fixed T and abundance fractions). Since 
halos with virial temperature T coo i must satisfy i coo ; ~ tu , 
see (5), T coo i has to get higher with time to increase 
the rate kij (which usually contains factors of the form 
exp(— Tij/T)). The sudden increase of T coo i at z ~ 30 
(log(l + z) ~ 1.5) is due to the decline of the fraction 
of molecular hydrogen which starts being destroyed by 
the radiation emitted by stars and quasars. As a conse- 
quence the main cooling process becomes collisional ex- 
citation cooling instead of molecular cooling. Since the 
former is only active at high temperatures (the coefficient 
rate k(T) contains a term exp(— 118348i£'/T') instead of 
exp(— 512K/T) for molecular cooling) the cooling tem- 
perature T coo i has to increase up to T coo i ~ 10 K. By 
definition T coo i is larger than the IGM temperature and 
usually much higher as can be seen in Fig. 3. However, at 
z ~ 9 when Tigm ~ 10 4 K is quite high due to reheating 
by the background radiation field we have T coo \ = Tigm 
since the IGM temperature is large enough to allow for 
efficient cooling. Then all virialized bound objects, with 
T > Tjqm, form baryonic clumps which can cool. The 
temperature T m represents a mass-averaged temperature: 
the matter within the IGM is associated to Tigm while 
virialized objects (hence with T > Tigm) are character- 
ized by a temperature defined as Min(T, 10 6 K). Since T m 
does not enter any of our calculations used to obtain the 
redshift evolution of the universe this crude definition is 
sufficient for our purpose which is merely to illustrate the 
difference between volume (Tigm) and mass averages. As 



can be seen from Fig. 3 we always have T m > Tigm as it 
should be. At large redshifts T m ~ Tigm since most of the 
matter is within the uniform IGM component, whereas at 
low redshifts z < 5 (log(l + z) < 0.8) the IGM tempera- 
ture declines as we explained previously while T m remains 
large since most of the matter is now embedded within col- 
lapsed objects where shock heating is important (and they 
do not experience adiabatic cooling due to the expansion). 

We show in Fig. 4 the cooling and heating times asso- 
ciated with various processes for the IGM as well as for 
the smallest halos M coo i which can cool at z. 
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0.5 1 1.5 

lttg(l+z) 




0.5 1 1.5 2 

log(l+z) 



Fig. 4. The cooling and heating times associated with the var- 
ious relevant processes in units of s tii(z) for the IGM (upper 
figure) and the halos defined by T coo i (lower figure). The labels 
are as follows: 1) collisional excitation, 2) collisional ionization, 
3) recombination, 4) molecular hydrogen, 5) bremsstrahlung, 
6) Compton, 7) photoionization heating and 8) cooling due to 
expansion (only for the IGM, see text). 

We can see in the upper panel that for large and 
small redshifts, z > 24 (log(l + z) > 1.4) and z < 5 
(log(l + z) < 0.8), all time-scales associated with the 
IGM are larger than the Hubble-time which means that 
the IGM temperature declines due to the adiabatic cool- 
ing entailed by the expansion of the universe. However, 
at intermediate redshifts z ~ 18 (log(l + z) ~ 1.3) the 
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smallest time-scale corresponds to heating by the back- 
ground radiation (theat) which means that Tigm increases 
during this period. Next, at z ~ 9, the IGM temper- 
ature becomes large enough to activate collisional exci- 
tation cooling so as to reach a temporary equilibrium 
where t coo i ~ theat while Tigm remains constant. Then, 
as we shall see later the universe gets suddenly reionized 
at z ri = 6.8 (log(l + z) = 0.9). This means that theat 
increases sharply as uhi declines (as well as nu e i and 
nHeii) as can be seen from (29). The cooling time due 
to collisional excitation follows this rise as the medium 
remains in quasi-equilibrium while the temperature de- 
clines slightly (the strong temperature-dependent factors 
like exp(— 118348i , T/T) in t coo i ensure it immediately ad- 
justs to theat, moreover these cooling rates are also pro- 
portional to Hh i, nHei and njieii) until both heating and 
cooling time-scales become larger than the Hubble time. 
Then this quasi-equilibrium regime stops as the medium 
merely cools because of the expansion of the universe. 
These various phases, which appear very clearly in Fig. 4, 
explain the behaviour of Tigm shown in Fig. 3 which we 
described earlier. The peak at z ~ 19 (log(l + z) ~ 1.3) 
of tcompton in the upper panel (curve 6) corresponds to 
the time when its sign changes (hence t Compton = 0). At 
higher redshifts Tigm is lower than the CMB tempera- 
ture (due to adiabatic cooling by the expansion of the 
universe), so that the gas is heated by the CMB photons 
while at lower z the IGM temperature is larger than Tcmb 
(due to reheating) so that the gas is cooled by the inter- 
action with the CMB radiation. 

The lower panel shows the cooling and heating times 
associated with the halos T coo i- As we have already ex- 
plained we can see that at high redshifts the main cooling 
process is molecular hydrogen cooling. Note however that 
for the IGM this process is always irrelevant. This dif- 
ference comes from the fact that the larger density and 
temperature of these virialized halos allow them to form 
more molecular hydrogen than is present in the IGM, so 
that molecular cooling becomes efficient. This was also 
described in detail in Tegmark et al.(1997) for instance. 
Of course at these redshifts we have t coo i,mol — s iff by 
definition of T cool . Then, at z < 27 (log(l + z) < 1.4) as 
molecular hydrogen starts being destroyed by the back- 
ground radiation the main cooling process becomes col- 
lisional excitation. Note that the corresponding cooling 
time gets smaller than s ta because the medium is also 
heated by the radiation so that the actual cooling time re- 
sults from a slight imbalance between cooling and heating 
processes. The sharp decrease of the various time-scales 
at z ~ 9 corresponds to a sudden increase of T coo i due 
to the rise of Tigm (which influences the cooling halos 
since T coo i > Tigm) also seen in Fig. 3 and in the upper 
panel of Fig. 4. Around z ^ 8 (log(l + z) ~ 0.9) we have 
tcooi *C theat and t coo i -C tn so that all virialized ha- 
los above Tigm can cool (T coo / = Tigm)- The feature at 
z ~ 5.3 (log(l + z) r-j 0.8) is due to reionization. 



Finally, we present in Fig. 5 the characteristic masses 
we encounter. The mass Md is obtained from (22): 

Air 

M d = (1 + A) igm PY R d ( 46 ) 

while M coo i corresponds to the halos which can cool at 
redshift z, as we have already explained. The mass M vir 
which follows closely the behaviour of Md describes the 
smallest virialized objects. It differs from Md because the 
density contrast is now A c instead of Aigm- The mass 
Mnl corresponds to the first non- linear scale defined by 
£ = 1. Note that after reionization Md ~ 3 10 6 Af Q while 
the usual Jeans mass would be Mj ~ lO 8 Af0. This is 
mainly due to the low density (1 + A)igm of the IGM, see 
(46) and Fig. 13 below. 
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Fig. 5. The redshift evolution of the characteristic masses Md, 
M V i r , M coo i and M NL in M©. 

By definition we have M d < M V i r < M CO ol- At z ~ 9 
we have M coo i = M V i r since T coo i — Tigm as we noticed 
earlier in Fig. 3. We also note that at large z our calculation 
is not entirely correct since our multiplicity functions are 
valid in the non-linear regime, for masses M <C M^l- 
However, at these early times the universe is nearly exactly 
uniform (by definition !) so that this is not a very serious 
problem. We can see that the first cooled objects which 
form in significant numbers are halos of dark-matter mass 
M ~ 1O 5 M which appear at z ~ 49 (log(l + z) - 1.7), 
when M coo i becomes smaller than Mnl. However, they 
only influence the IGM after z < 24 (log(l + z) < 1.4) 
when reheating begins. 

8.3. Reionization of the IGM 

After the radiation emitted by quasars and stars reheats 
the universe, as described in the previous section, it will 
eventually reionize the IGM. We present in Fig. 6 the evo- 
lution with redshift of the background radiation field and 
of the comoving stellar formation rate. Within the frame- 
work of our model the latter is a good measure of the 
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radiative output from galaxies, see (12), as well as from 
quasars, see (17), since we note that the quasar mass hap- 
pens to be roughly proportional to the stellar mass. 




0.5 1 1.5 2 

log(l + z) 




0.5 1 1.5 2 

log{l+z) 



Fig. 6. The redshift evolution of the UV flux J21 (upper panel) 
and of the comoving star formation rate dp a /dt (lower panel). 
The data points are from Giallongo et al.(1996) (square), 
Cooke et al.(1997) (filled square), Vogel et al.(1995) (trian- 
gle, upper limit), Donahue et al.(1995) (filled triangle, upper 
limit) and Kulkarni & Fall(1993) (circle). The dashed line in 
the lower panel shows the effect of the absorption of high en- 
ergy photons by the neutral hydrogen present in the IGM and 
in Lyman-a clouds. 

The upper panel of Fig. 6 shows the UV flux J21 in 
units of 10 -21 erg cm -2 s _1 Hz -1 sr _1 defined by (38). 
We can see that the UV flux rises very sharply at z ~ 6.8 
(log(l + z) ~ 0.9) which corresponds to the reionization 
redshift z r % when the universe suddenly becomes optically 
thin, so that the radiation emitted by stars and quasars 
at large frequencies is no longer absorbed and contributes 
directly to J u . This appears clearly from the lower panel. 
Here the solid line shows the comoving star formation rate, 
obtained from (10), while the dashed line shows the same 
quantity multiplied by a luminosity-weighted opacity fac- 
tor exp(— tl) which describes the opacity due to the IGM 
and Lyman-a clouds (see below (51), (52) and Fig. 10). 



Thus, we can see that while the star formation rate evolves 
rather slowly with z the absorption term varies sharply 
around z r i- Hence the universe is suddenly reionized at 
z r i (when the ionized bubbles overlap: Qhii — 1) on a 
time-scale very short as compared to the Hubble time be- 
cause of the strongly non-linear effect of the opacity. We 
note that for z < 1 our star-formation model is some- 
what simplified, as we explained earlier, because we de- 
fined all galaxies by a constant density contrast A c while 
cooling constraints should be taken into account, as in VS 
II, and we also used approximations for the stellar content 
of galaxies which are not strictly valid for all galactic halos 
at these low redshifts. The reader is referred to VS II for a 
more precise description of the low z behaviour. However, 
our present treatment is sufficient for our purposes and 
still provides a reasonable approximation at z < 1. 

We display in Fig. 7 the background radiation spectrum 
J„ at four redshifts: z\ = 7.3 (before reionization), Z2 = 
6.4 (after reionization), 2:3 — 3 and z 4 = 0. 




log(f) [ eV ] 

Fig. 7. The background radiation spectrum J„ (in units of 
10~ 21 erg cm -2 s _1 Hz -1 sr _1 ) at the redshifts z\ = 7.3 (solid 
line, prior to reionization), 22 = 6.4 (lower dashed line, after 
reionization), Z3 = 3 (upper dashed line) and 24 = (solid 
line) . 

We see that the ionization edges corresponding to HI, 
Hel and Hell can be clearly seen at high z before the uni- 
verse is reionized. Then the background radiation is very 
strongly suppressed for v > 13.6 eV due to HI and Hel 
absorption. Of course at very large frequencies v il 1 keV 
where the cross-section gets small we recover the slope 
J v oc v~ a of the radiation emitted by quasars. At low red- 
shifts z < z r i after reionization the drop corresponding to 
Hel disappears as Hel is fully ionized and its number den- 
sity gets extremely small, as we shall see below in Fig. 9. 
However, even at z ~ 3 the ionization edges due to HI 
and Hell are clearly apparent and J„ is still significantly 
different from a simple power-law. At low redshifts z ~ 
the background radiation is much smoother since its main 
contribution comes from radiation emitted while the uni- 
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verse was ionized and optically thin. However, its inten- 
sity is smaller than at z ~ 4 because the quasar luminosity 
function drops at low z, see Fig.l, while the universe keeps 
expanding. 

We show in Fig. 8 the redshift evolution of the ioniza- 
tion and recombination times ti on and t rec of the IGM, 
divided by tn ■ 
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Fig. 8. The redshift evolution of the ionization and recombi- 
nation times ti 0n (solid line) and t rec ,iGM (dashed line) of the 
IGM, divided by the Hubble time ts- The horizontal solid line 
only shows tn for reference. The recombination times t T ec,mean 
(uniform medium) and t rec ,bubbie (ionized bubbles) are defined 
in the main text. 



More precisely, the ionization time ti on is defined by: 
C = J ^ ^ °HiV) dv (47) 



while the recombination time within the IGM is: 



t rec,IGM = Ot{T IGM ) C n U e 



(48) 



where a is the recombination rate, C n the clumping factor 
and n e - the mean electron number density, from (27) and 
(26). We also display for reference the recombination time 
which would correspond to a uniform medium with the 
mean density of the universe: 



t 



rec.mean 



(1 + A) n C n t rec JGM 



(49) 



Finally, we show the recombination time within ionized 
bubbles (where all hydrogen atoms are ionized): 



t 



rec.bubble 



a (3 10 4 K) C n n H 



(50) 



The recombination time grows with time at high z as the 
mean density decreases with the expansion, although this 
is somewhat balanced by the increase of the clumping fac- 
tor C n (see below Fig. 13). In particular, the decrease of 
tree, igm around z ~ 30 (log(l + z) ~ 1.5) is due to the 
growth of C n . The sharp drop at z ~ 7 (log(l+z) ~ 0.9) is 



due to reionization which suddenly increases the number 
density of free electrons. After reionization the recombina- 
tion time characteristic of the IGM keeps decreasing (while 
t r ec,mean increases slightly since the density declines) be- 
cause of the growth of the clumping factor C n , see Fig. 13, 
which overides the decline of the mean universe density. 
The recombination time within ionized bubbles t reCjU ubbie 
follows the change of the mean universe density and of the 
clumping factor C n ■ At large z it is much smaller than the 
mean IGM recombination time since the IGM is close to 
neutral. At low z it becomes larger than t re c,iGM since 
the IGM is suddenly reionized with a temperature Tiqm 
which declines after reheating and gets lower than 3 10 4 K, 
see Fig. 3. The ionization time is very large at high z since 
the UV background radiation is small. Then it decreases 
very sharply at z r i when the universe is reionized and the 
background radiation suddenly grows as the medium be- 
comes optically thin, as seen in Fig. 6. The reionization red- 
shift corresponds to the time when ti on becomes smaller 
than tfj, somewhat after it gets smaller than t reCl i G M- 
Thus t rec jgm does not play a decisive role since it is never 
the smallest time-scale around reionization. 

Finally, we show in Fig. 9 how the chemistry of the IGM 
evolves with time as the temperature Tjqm and the UV 
flux J21 vary with z. 

We see very clearly in the upper panel the redshift of 
reionization z„ = 6.8 (log(l + z) = 0.9) when the fraction 
of neutral hydrogen nni jnu declines very sharply while 
the UV flux J21 suddenly rises, as was shown in Fig. 6. We 
note that at low redshifts z < z r i the neutral hydrogen 
fraction decreases more slowly down to rim fun ~ 10 -8 . 
The fractions of electrons and ionized hydrogen start in- 
creasing earlier at z ~ 19 (log(l + z) ~ 1.3) but of course 
they remain small until z ri . The abundance of molecu- 
lar hydrogen decreases sharply at a rather high redshift 
z ~ 27 (log(l + z) oj 1.4) due to the background radia- 
tion, as we noticed on Fig. 4. The lower panel shows that 
helium gets fully ionized simultaneously with hydrogen. 
In particular, although there remains a small fraction of 
Hell (~ 10~ 4 ) the abundance of Hel gets extremely small. 
We note that at low redshifts z < 2 the fraction of Hell 
does not evolve much (and even slightly increases) while 
the HI abundance keeps declining. This is due to the fact 
that the radiation relevant for helium ionization comes 
from quasars whose luminosity function drops at low z 
as shown in Fig.l while an important contribution to the 
hydrogen ionizing radiation is provided by stars and the 
galaxy luminosity function declines more slowly with time 
at low z, as seen in Fig. 6 or in VS II. We shall come back 
to this point in Sect. 8. 7. 

8.4- Opacities 

As we explained previously the radiation emitted by stars 
and quasars at high frequencies (v > 13.6 eV) is absorbed 
by the IGM and discrete clouds as it propagates into the 
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Fig. 9. The redshift evolution of the chemistry of the IGM. 
The upper panel shows the ionization state of hydrogen, as 
well as the fraction of molecular hydrogen and electrons. The 
lower panel presents the ionization of helium. 

IGM. This leads to the extinction factors exp(— r s (x)) and 
exp(— tq(x)) in the evaluation of the source terms (13) 
and (21) for the background radiation. We define here the 
"luminosity averages" t\ gm nhi for both continuous and 
discrete components by: 

•-^-(tr/*Mf (^)w^ w (ei) 

where t s (x) is the relevant opacity (from the IGM or 
clouds for sources x, see (44) ) at the frequency 20 eV 
below the Hel ionization threshold without taking into 
account the factors Qui which describe ionized bubbles. 
The subscript s refers to the fact that we consider here 
the opacity which enters into the calculation of the stellar 
radiative output. The quasar-related opacity mainly dif- 
fers through the factor Aq, see (42). The weight dM s /dt 
corresponds to a luminosity weight as we noticed earlier, 
see (12), (17) and (6). 

We show in the upper panel of Fig. 10 the redshift evo- 
lution of the opacity from the IGM (tj GM , dashed line) 
and from "Lyman-a clouds" (t^ hi , solid line). We can see 
that both contributions have roughly the same magnitude 



Fig. 10. Upper panel: the redshift evolution of the opacity 
from the IGM (dashed line) and "Lyman-a clouds" (solid line) 
which enters the absorption factors in the calculation of the 
radiative output from stars and quasars. Lower panel: evolution 
of the filling factors Qhii q (ionized bubbles around quasars, 
solid line), Qhii b (around galaxies, upper dashed line) and 
Qmi.cioud (lower dashed line). 

before reionization, except at very high redshifts z > 50 
(log(l + z) > 1.7) when very few structures exist as shown 
in Fig. 12. Prior to reionization the opacity is large and the 
background radiation quite small, as seen in Fig. 6. At z r i 
the opacity suddenly declines while J21 rises sharply, due 
to the strong non- linear coupling between r and J v , as we 
explained in the lower panel of Fig. 6 where we presented 
the influence of the total opacity tl- 

T~L = T IGM + T NHI (52) 

At low redshifts z < z r i when the universe is reionized the 
opacity due to discrete clouds becomes much larger than 
the IGM contribution (although it is very small) because 
the density of neutral hydrogen is now proportional to the 
square of the baryonic density (in photoionized regions) 
and most of the matter is embedded within collapsed ob- 
jects. 

The opacities t l were shown in the upper panel of 
Fig. 10 without the filling factors Qhi which enter the 
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actual evaluation of the source terms (13) and (21), see 
(40). The filling factors Qhii, describing the volume frac- 
tion occupied by ionized bubbles, are shown in the lower 
panel of Fig. 10. We can see that Qhii 3 and Qhii q increase 
with time as structures form and emit radiation while the 
IGM density declines. When these ionized bubbles over- 
lap (Qhii = 1) the universe is reionized. At low redshifts 
the coefficient Qhii, cloud declines because the background 
radiation is large while the quasar number density drops 
(indeed Qhii. cloud measures the volume occupied by the 
"spheres of influence" of quasars). 

Next, we can evaluate the mean opacities thi{z) and 
THeii(z) seen on a random line of sight from z = to 
a quasar located at redshift z. We present in Fig. 11 the 
contributions from both the uniform IGM component and 
the discrete Lyman-a clouds. 
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Fig. 11. The redshift evolution of the average opacities thi 
and t Heii along a random line of sight produced by "Lyman-a 
clouds" . The dashed lines show the opacities from the uniform 
IGM. The data points are from Press et al.(1993) (circles), Zuo 
& Lu (1993) (filled circles) for hydrogen, and from Davidsen et 
al.(1996) (filled rectangle) and Hogan et al.(1997) (cross) for 
helium. 

At high redshifts prior to reionization the main contri- 
bution to the opacity is provided by the IGM which con- 
tains most of the matter. However, at low z when the IGM 
is ionized most of the absorption comes from the Lyman-a 
clouds. At large z the Hell opacity is very small because 
most of the helium is in its neutral form Hel. We refer the 
reader to Valageas et al. (1999a) for a much more detailed 
description of the properties of the Lyman-a clouds at 
low z. We can see that our predictions show a reasonable 
agreement with observations for the hydrogen opacity. At 
low redshifts z < 1 the influence of star-formation which 
consumes and may eject some of the gas (which we did 
not take into account here) could explain the relatively 
high opacity we obtain. The helium opacity we find is also 
close to observations. This is due to the fact that the UV 
radiation spectrum shows strong ionization edges, even at 



relatively low redshifts z ~ 3 (log(l + z) ~ 0.6), see Fig. 7. 
Hence the ratio NhbIi/Nhi is rather large, see Fig. 9, 
which explains why we get a better agreement than Zheng 
et al.(1998) for instance (see also Valageas et al. 1999a). 
In particular, we have nHdi /njjein ^ nni/nnn since 
n>Heii /nHein ~ 10 -4 while njji/nHii ~ 10 -8 at low red- 
shift. We note that the observed Hell opacity strongly con- 
strains the quasar contribution to the reionization process 
since stellar radiation is small at high frequencies (due 
to the near blackbody behaviour of stellar spectra). In 
particular, it implies that one needs a population of faint 
QSOs (Mb > —26.7) in order to reionize helium but z r i 
should not be too large so that there is still an apprecia- 
ble density of Hell. In other words, as we noticed above, 
the UV radiation field must still display strong ionization 
edges, which means that it has not had enough time to be 
smoothed out by the radiation emitted since z r i when the 
medium is optically thin. 

8. 5. Stellar properties 

Our model also allows us to obtain the fraction of matter 
within virialized or cooled halos, as well as in stars. 
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Fig. 12. The redshift evolution of the fraction of matter en- 
closed within virialized halos (F V i r ), cooled objects (F coo i) and 
stars (Fstar)- The lower solid line is the volume fraction F V i r<vo i 
occupied by virialized objects. 

We show in Fig. 12 the fraction of matter within virial- 
ized halos (Evir, upper dashed line), cooled objects (F coo i, 
upper solid line) and stars (F s t a r, lower dashed line). The 
first two quantities are simply obtained from (3). Of course 
we have: F star < F coo i < F V i r . Around z ~ 9 we note 
that F coo i = F vir because as we explained previously at 
this time all virialized objects (with T > Tjqm) can cool 
efficiently (T coo i = 7/Gm)- The fraction of matter within 
virialized halos increases very fast at high redshifts z ~ 49 
(log(l + z) ~ 1.6) as M vir becomes smaller than M^^, 
see Fig. 5, when dark matter structures form on scale Rd- 
However, until z ~ 15 (log(l + z) ~ 1.2) the mass within 
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cooled halos remains much smaller because cooling is not 
very efficient so that T coo ; ^> Tiqm, see Fig. 3. At low red- 
shifts z < 5 (log(l + z) < 0.8) both F vir and F coo i get 
close to unity since most of the matter is now embedded 
within collapsed and cooled halos (even though T coo i be- 
comes again much larger than Tigm- we are so far within 
the non-linear regime that even T coo i is small compared 
to the characteristic virial temperature of the structures 
built on scale Rd). Of course the mass within stars grows 
with time, closely following F coo i. Note however that it 
is not strictly proportional to F coo i since an increasingly 
large fraction of the gas within galaxies is consumed into 
stars. The fraction of volume F V i r ^ vo i occupied by virial- 
ized objects always remains small as it satisfies: 
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Fig. 13. Upper panel: the redshift evolution of the clump- 
ing factors Ct and C„- Lower panel: the overdensities 
(1 + A)igm(z) characteristic of most of the volume of the uni- 
verse at redshift z at scale Rd and (1 + A) n . 

We show in the upper panel of Fig. 13 the clumping fac- 
tor Cb defined by (25). The expression (25) shows clearly 
that at large redshifts where there are very few collapsed 
baryonic structures C'b — 1 while at low z when most 
of the baryonic matter is within virialized halos (note 
this is always true for dark matter on sufficiently small 
scales) we have ~ A c . We can see in the figure that 



Cb usually increases with time as the hierarchical cluster- 
ing process goes on. The temporary decrease at z ~ 15 
(log(l + z) ~ 1.2), which also appears in the lower panel 
and in Fig. 12, is due to the reheating of the universe, 
shown in Fig. 3, which increases the "damping" length Rd 
and mass scale Md, as seen in Fig. 5. As a consequence, 
small objects which were previously well-defined entities 
suddenly see the IGM temperature become larger than 
their virial temperature. Hence they cannot retain effi- 
ciently their gas content and they lose their identity. We 
note that neglecting the clumping of the gas would lead 
to a higher reionization redshift z r i — 7.3 since it would 
underestimate the efficiency of recombination. The clump- 
ing factor C n displays a behaviour similar to Cf, but it is 
usually smaller since it does not include the deep halos 
which can cool. We display in the lower panel of Fig. 13 
the overdensity (1 + A)igm(z) characteristic of most of 
the volume of the universe at redshift z when seen on scale 
Rd- While structures form and the matter gets embedded 
within overdensities which occupy a decreasing fraction of 
the volume (when seen at this scale Rd) the "overdensity" 
(1 + A)iqm which characterises the medium in-between 
these objects (halos or filaments) declines. The density 
contrast (1 + A)„ which corresponds to the IGM and shal- 
low potential wells which do not form stars (but constitute 
Lyman-a clouds) decreases more slowly since it only ex- 
cludes the high virial temperature halos with T > T coo /. 




log(H-z) 

Fig. 14. The redshift evolution of the metallicities Z c 
(star- forming gas), Z s (stars), Zh (galactic halos) and Z m 
(matter average). The data points are from Pettini et al.(1997) 
for the zinc metallicity of damped Lyman-a systems. 



We present in Fig. 14 the redshift evolution of the mean 
metallicities (in units of solar metallicity) Z c (within the 
star-forming gas located in the inner parts of galaxies, 
upper solid line), Z s (within stars, upper dashed line) and 
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Zh (within galactic halos, lower dashed line). We use the 
mass average over the various galactic halos: 



JZ(M) M 9 (M)^ 



(54) 



where fi g (M)dM/M is the galaxy mass function defined 
as in (3) . The reader is referred to VS II for a detailed de- 
scription of these metallicities (note that we only consider 
here the abundance of Oxygen or any other element that is 
mainly produced by SN II since we did not include SN I in 
our model). The lower solid line corresponds to a "matter 
averaged" metallicity Z m defined by Z m = F coo iZh- Thus, 
although we do not include explicitly in our model any 
contamination of the IGM by heavy elements produced 
within galaxies, Z m defined in this way provides an up- 
per bound for the mean IGM metallicity (corresponding 
to very efficient mixing). If galaxies do not eject metals 
very deeply within the IGM its metallicity could be much 
smaller. The mean metallicity of Lyman-a clouds associ- 
ated to galactic halos (limit or damped systems) is Zh- Our 
results agree well with observations by Pettini et al.(1997) 
for damped Lyman-a systems. Note that there is in fact a 
non-negligible spread in metallicity over the various halos. 

8. 6. Consequence for the CMB radiation 

After reionization, CMB photons may be scattered by 
electrons present in the gas. We write the corresponding 
Thomson opacity up to a redshift z as: 



r dt 



a T n e (z) 



(55) 



where n e (z) is the mean electron number density at red- 
shift z. We take: 



n e (z) = (l-Y) 



n b p(z) 



Tie 

n H J IGM 



(56) 



which means that we use the same electron fraction in 
clouds as for the IGM (note that we calculate the IGM 
electron number density together with the ionization state 
of hydrogen and helium). Then, CMB anisotropies are 
damped on angular scales smaller than the angle sub- 
tended by the horizon at reionization. We use the analytic 
fit given by Hu & White (1997) to obtain the damping fac- 
tor Rf of the CMB power-spectrum Ci from the optical 
depth T es . The results are shown in Fig. 15. 

We can check that the total opacity T es ~ 0.023 is quite 
small because reionization occurs rather late at z r i ~ 6.8. 
This also implies that the damping factor Rf remains 
close to unity: Rf ~ 0.95 for large I. Another distor- 
tion of the CMB radiation is the Sunyaev-Zeldovich effect 
which transfers photons from the Rayleigh-Jeans part of 
the spectrum to the Wien tail when they are scattered 




Fig. 15. Upper panel: the optical depth r es for electron 
scattering. Lower panel: damping factor Rf for the CMB 
power-spectrum. 



by hot electrons. The magnitude of this perturbation is 
conveniently described by the Compton parameter y: 



y{z) 



-f 

Jo 



dt , kT 

c—dz ar n e t 

dz m e c z 



(57) 



We can first consider the contribution of the IGM gas, 
using in (57) the temperature and the electronic density of 
this uniform component. Then, we estimate the distortion 
due to the hot gas embedded within virialized objects. We 
can write this latter contribution as: 



Vhalos — / dyiGM 



1 



T 

J. 77- 



(1 + A)/ GM Tjgm 



(58) 



where we used the same electronic fraction for halos and 
the IGM. The temperature T mvir is the "mass averaged" 
temperature of virialized objects while F vir is the mass 
fraction within collapsed objects displayed in Fig. 12. The 
results are presented in Fig. 16. 

The Compton parameter yiGM due to the IGM first in- 
creases rather fast with z until reionization, together with 
n e and Tjgm- After z r i it reaches a plateau and does not 
grow any more since at these large redshifts the universe 
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Fig. 16. The Compton parameter y up to redshift z describing 
the Sunyaev-Zeldovich effect from the IGM (dashed line) and 
virialized halos (solid line). 



is almost exactly neutral. The contribution ytaios of viri- 
alized objects is much larger at low z than yiGM since the 
temperature of these collapsed halos is much higher than 
Tigm, as shown in Fig. 3. We can note however that yuaios 
grows much slower and becomes close to its asymptotic 
value earlier than yiGM- This is due to the fact that the 
characteristic temperature of virialized halos declines at 
larger z, see Fig. 3, and the mass fraction they contain also 
decreases (while the IGM undergoes the opposite trends). 
A more detailed description of the Sunyaev-Zeldovich ef- 
fect due to clusters, and its fluctuations (which in fact 
have the same magnitude as the mean), will be presented 
in a future article. 

8. 7. Contributions of quasars and stars 

In our model the radiation which reheats and reionizes 
the universe comes from both quasars and stars. At large 
frequencies v > 24.6 eV most of the UV flux is emitted 
by quasars so that stars play no role in the helium ion- 
ization. However, at lower frequencies both contributions 
have roughly the same magnitude. We present in Fig. 17 
the redshift evolution of the radiative output due to stars 
and quasars. 

Thus, we dehne the "instantaneous" radiation fields: 



= ±t H S u 



(59) 



where tjj(z) is the Hubble time, see (30). From these quan- 
tities we define the averages Jhi-, s ,q( z ) an d Jhei:s,Q( z ) 
as in (38). This provides a measure of the radiative out- 
put above the ionization thresholds vjji and vjjei due 
to stars and quasars. We can also derive the ionization 
times t\ s q as in (47). We can see in the upper panel 
that at reionization z ~ z, ri the contributions to HI ion- 
izing radiation from stars and quasars are of the same 
order. However, since the quasar spectrum is much harder 
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Fig. 17. Upper panel: the redshift evolution of the "instan- 
taneous" UV fluxes J l HI and J}f e i due to stars (solid lines) 
and quasars (dashed lines). The dotted line is the UV flux 
J21 shown in Fig. 6. Lower panel: the "instantaneous" ioniza- 
tion times t\ on . s q due to stars (solid line) and quasars (dashed 
line). The dotted curve is the recombination time as in Fig. 8 
while the horizontal solid line is the Hubble time tn- 



than stellar radiation we have J\j e j.Q 3> J 



Hel:. 



so that 



quasars are slightly more efficient at reheating and reion- 
izing the universe (the additional factor (v — Vj) in (29) 
increases the weight of high energy photons which also 
remain longer above the threshold vhi while being red- 
shifted). At low z we can see that the quasar radiative 
output declines much faster than the stellar source term. 
Of course, this is due to the sharp drop at low redshifts 
of the quasar luminosity function. This decrease of S u q 
as compared to S us comes from two effects: i) as time 
increases the "creation time-scale" of halos of the rele- 
vant masses M ~ 10 12 M Q (through merging of smaller 
sub-units, measured by tg/tM) grows and ii) there is less 
gas available to fuel the quasars (which even disappear) 
while old stars can still provide a non-negligible luminosity 
source for galaxies. We can note in the upper panel that 
the redshift evolution of the background radiation field 
J2i{z) actually produced by stars and quasars does not 
exactly follow the "instantaneous" quantities J\jj. s q{z) 
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since one must take into account the expansion of the uni- 
verse and deviations from equilibrium. This explains the 
slower increase of J21 at well as the relatively 

low approximate ionization times t l ion . s q at this epoch. 

9. Critical universe 

We now consider the case of a critical universe Q = 1 
with a CDM power-spectrum (Davis et al.1985) normal- 
ized to er 8 = 0.5. We also choose = 0.04 and H = 60 
km/s/Mpc. Thus, as in the previous case of an open uni- 
verse our model is consistent with the studies described in 
VS II and Valageas et al.(1999a). 

9.1. Quasar luminosity junction 

We first present our results for the redshift evolution of 
the quasar luminosity function in Fig. 18. 

We can check that our results are similar to those ob- 
tained previously for an open universe and they again 
agree reasonably with observations. This is not very sur- 
prising since we use the same physical model so that we re- 
cover a similar behaviour. We present in Fig. 19 the quasar 
number counts we obtain from our model. 

We can see that our results are similar to those dis- 
played previously in Fig. 2 and that our predictions are 
still marginally consistent with the lack of observation in 
the HDF. 

9.2. Reheating 

We present in Fig. 20 the redshift evolution of the IGM 
temperature Tjgm, the virial temperature T coo i of the 
smallest objects which can cool at a given time, and the 
mass averaged temperature T m . 

We can see that our results are again very close to 
those obtained for an open universe. Indeed, the struc- 
ture formation process is quite similar and it must agree 
with the same observations (quasar and galaxy luminosity 
functions, Lyman-a column density distribution) at low z. 




11 12 13 14 15 



log(L B /L s ) 

Fig. 18. The evolution with redshift of the B-band quasar 
luminosity function in comoving Mpc -3 , as in Fig.l. The data 
points are from Pei (1995). 



9.3. Reionization 

We display in Fig.21 the redshift evolution of the back- uncertainties involved in the galaxy and quasar formation 

ground radiation and the comoving star formation rate. Presses are probably too large to favour significantly 

We can check that we recover the behaviour obtained one of these tw0 P ossible scenarios ( as compared to the 

previously for an open universe. However, the reionization other )' However, both models are consistent with present 

redshift z n = 5.6 is smaller than previously. This is related observations, 
to the lower normalization as of the power-spectrum as 

compared to the previous case. This leads to fewer bright -^n r; onc i us i on 
quasars at high z (compare Fig. 18 and Fig.l) and to a 

smaller radiative output. In this article we have described an analytic model for 

We can also check that the hydrogen and helium reion- structure formation in the universe which deals simultane- 

ization process is close to our previous results (as for the ously with quasars, galaxies and Lyman-a clouds, within 

reheating). Thus, for most practical purposes both criti- the framework of a hierarchical scenario. This allows us 

cal and open cosmologies allow reasonable reheating and to study the reheating and reionization history of the 

reionization histories which are very similar. In fact, the universe consistently with the properties of these vari- 
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Fig. 19. The quasar cumulative V-band number counts. The 
dashed line shows the counts of quasars with a magnitude 
brighter than V located at redshifts 3.5 < z < 4.5 while the 
solid line corresponds to 3.5 < z < z ri . 
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Fig. 20. The redshift evolution of the IGM temperature Tigm 
(solid curve), the virial temperature T coo i (upper dashed curve) 
and the mass averaged temperature T m (lower dashed curve), 
as in Fig. 3. 

ous classes of objects. We have shown that for both a 
critical and an open universe our predictions agree rea- 
sonably well with observations. However, as was noticed 
by Haiman et al.(1998) it appears that the observational 
constraints on the quasar luminosity function are already 
strong. Moreover, the Gunn-Peterson test for Hell pro- 
vides stringent additional constraints on the quasar con- 
tribution to the UV radiation field and on the reionization 
redshift. Thus, although our model in its simplest version 
(i.e. as described here, with no additional cutoffs for the 
quasar multiplicity function) is marginally consistent with 
the data, further observations of the helium opacity and of 
the quasar number counts (e.g. with the NGST) could pro- 
vide tight constraints on such models where reionization 
is produced by QSOs. On the other hand, since reion- 
ization occurs rather late z rl < 7 the damping of CMB 
anisotropics is quite small. 
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Fig. 21. The redshift evolution of the UV flux J21 (upper 
panel) and of the comoving star formation rate dp„/dt (lower 
panel) for the case of an open universe. The dashed line in the 
lower panel shows the effect of the absorption of high energy 
photons by the neutral hydrogen present in the IGM and in 
Lyamn-a clouds. 

We can note that our predictions are similar to some 
results obtained by Gnedin & Ostriker (1997) with a nu- 
merical simulation (but for a different cosmology). More- 
over, the fact that our model agrees reasonably with ob- 
servations for Lyman-a clouds, galaxies, quasars and con- 
straints on the reionization process, strongly suggests that 
its main characteristics are fairly realistic. Thus, it pro- 
vides a simple description of structure formation in the 
universe, from high redshifts after recombination down to 
the present epoch. 
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